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ABSTRACT 

A general approach is documented as a guide to aid in the formulation 
and implementation of on-line, real time computer simulations* A computer 
program MULNUCl, is developed as an on-line, real time computer simu- 
lation of antisubmarine warfare in a multiple burst nuclear environment. 
The principals of the game are a submarine armed with torpedoes, and two 
destroyers equipped with stand-off antisubmarine weapons. The simu- 
lation is intended as a demonstration of the on-line capabilities of the 
United States Naval Postgraduate School computer system and as a tool 
for further study of the factors involved in a representative ASW 
operational environment. 



9 



ua 

TABLE OF CONTENTS 

Section 

1. Introduction 

2. Formulation of On-Line Simulations 

Initializing 
An Iterative Loop 
A critique 

3. Program Modularization 

Advantages 
F lexibility 

4. System Variables, Coordinates and Notation 

5. MULNUCl - An Example 

Executive Control 
Plot Generator 
Interactions 
Sonar Contact Model 
Contact Tracking Model 
Weapon Firing Model 
Evaluation Model 
Radiation Model 
Submarine Logic Model 
Sequencing of Blocks 
Man Machine Interface 
Timing 

6. Conclusions and Acknowledgements 
Bibliography 



DUDLEY KNOX LiORARV 

NAVAL POSTGRADUAi E SC! 
MONTEREY CA 93943-5101 



Page 

7 

9 

9 

10 

11 

12 

12 

13 

14 
16 
18 
21 
24 

24 

25 

25 

26 

26 

27 

28 
30 

36 

38 

39 



3 



Section Page 

Appendix I - List of Variables 40 

Appendix II - Logic Flow Diagrams 50 

Appendix III - Subroutines and CDC 160 Executive Routine 69 

Appendix IV - Program Listing 74 



LIST OF ILLUSTRATIONS 



Figure 

1. Schematic Layout of Computer Complex 

2. Typical Left Tube Display on DD 65 

3. Typical Right Tube Display on DD 65 

4. DD 65 Console Arrangement 

5. Overlay used on DD 65 Keyboard 



Page 

17 

22 

23 

33 

35 



5 



1. INTRODUCTION 



Simulation is a useful tool of the operations analyst. This is 
not a new concept, the first recorded simulations were conducted by 
the Chinese in the form of "war chess" and were probably used to teach 
young men some of the concepts of battle without the inherent danger 
of loss of life. Later accounts of simulations have been recorded by 
the Prussians, French and Germans [7]. 

This technique of analyzing a problem by simulation is now employed 
by all branches of the scientific community. The present day high speed 
digital computer has given rise to a rapid expansion in the use of 
simulation as a method of solution to military, scientific, management 
and many other types of problems. 

Convincing the reader or observer that a simulation "models the 
real world" is one of the primary problems confronting the analyst 
who uses simulation techniques in the solution of a problem. One way 
to minimize this doubt is to increase the role of the human in the 
simulation. This can, in many cases, be done by the techniques of on- 
line simulation. While on the one hand, on-line simulation increases 
the complexity of the problem by nature of the man machine interface, 
at the same time, on-line simulation adds very complex logic (the man) 
to the problem with minimal effort on the part of the designer. 

The design of a simulation, and in particular an on-line simulation, 
can appear to be a formidable task when first considered. However, if 
the designer has an approach in mind and proceeds in an organized 
manner, the problem usually divides into small parts that are each 
relatively simple. It is the aim of this thesis to take a representative 
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problem and develop a computer simulation that can be used as a 



reference for the construction of on-line, real time, computer simu- 
lations in general. 
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2. FORMULATION OF ON-LINE SIMULATIONS 



In the formulation of on-line simulations the designer must first 
Lay down a good foundation in the form of a well-pianned outline* This 
outline must begin by initializing and setting the scene for the simu- 
lation. When this has been accomplished the designer must turn his 
attention to the development of a loop that will include all the actions 
and interactions expected to occur in the simulation. Included as an 
integral part of this loop is a timing mechanism that is flexible enough 
to allow the simulation to proceed at any rate required. A third part 
of the formulation is concerned with providing a critique of the simu- 
lation, either as a running critique or a compilation of pertinent facts 
at the end of the simulation. 

INITIALIZING 

Initializing is the term applied to that portion of the simulation 
which is executed before "play* 1 begins. It includes data input, assign- 
ment of particular values to the parameters and entering starting values 
required for indexing the logic. Herein is provided the flecibility 
required of any useful simulation. In the initializing portion the ground 
work must be formed for performing sensitivity analysis if such analysis 
is required. The initializing portion of the simulation must allow enough 
flexibility to provide for the various scenarios possible in the particular 
simulation. Clearly then, the initializing portion must be designed with 
these purposes as primary criteria and subject to boundary conditions, 
such as equipment capabilities. 
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AN ITERATIVE LOOP 



In general, many on-line simulations contain an iterative loop that 
is cycled for each time period. Therefore, one of the first considerations 
to be made in this phase of the design is the determination of the stepping 
interval. Several factors are involved in this selection, the most 
important being the assurance that the stepping interval is compatible 
with the logic flow of the situation being simulated. The designer must 
then consider the amount of time required to accomplish the most com- 
plicated simulation situation that can occur in one cycle. These con- 
siderations complete, the simulator must insure that the cycling is such 
that the player is not bored with the data/action as presented and also 
that this data/action is not presented at a rate too rapid for the player 
to fully comprehend. With these considerations in mind, a tentative 
looping cycle can be constructed and the designer may continue to develop 
the necessary routines required to complete the iterative loop. 

The loop should now be fashioned in its most elementary form. The 
designer must consider several tasks that must be accomplished during 
each cycle. These include, for example, advancing all participants one 
time cycle, consideration of the possible interactions that can occur due 
to these moves, tabulation of the results of such interactions, presen- 
tation of output to the player, permitting the player to communicate with 
the simulation, delaying the next cycle until the proper time interval 
has transpired, and possible other considerations dependent upon the 
particular simulation. 

Care and planning must be exercised in the construction of this loop 
since this is the foundation upon which the designer is to build his 
simulation. If logical errors appear in the order of these routines or 
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a component is not considered in this loop, the remaining portion of 
the design will be difficult, if not impossible. Planning in this 
portion of the development will be time well spent. 

THE CRITIQUE 

The critique is that portion of the simulation in which the entire 
simulation, or any integral part, is analyzed and the results compiled 
in condensed form. There are several basic tecnniques that may be 
explored in this part of the simulation. The least complex of which is, 
in most cases, a complete "recording" of the game that can be "replayed" 
at a later time at any speed desired. A more complex approach to the 
problem of performing a critique of a simulation is that of including 
a recording/analyzing routine in the iterative loop. This routine would 
extract the desired information during each cycle of the loop. It is 
apparent that not every cycle need be recorded, and therefore, a decision 
logic must be included that will extract all necessary information. 
Associated with this technique must be a recording routine that can be 
queried at intervals or at the conclusion of the simulation. 

Much care must be exercised in formulating the design and location 
of the critique routine since this is the major mode of analysis available 
to the analyst. Flexibility and adaptability are considerations that 
must be made to allow che simulation to be fully appreciated as an 
analytical tool and not just a "parlor game." 
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3. PROGRAM MODULARIZATION 



In general, programmers attempt to modularize their programs. There 
are several reasons for this, the most obvious being that of providing 
logical grouping of ideas. The large scale simulation is usually modu- 
larized by the use of subroutines. This use of subroutines is convenient 
because: 

a. Several separate groups may be working on various sections of 
the problem, and in many instances the use of subroutines is the best 
technique. 

b. Computations which are to be called upon several times in the 
main program are best handled by the use of subroutines. 

c. The program may be of a magnitude such that the entire program 
cannot be compiled in one pass. 

These reasons are valid for the large scale simulation. In simu- 
lations which are moderate to small in size the use of subroutines may 
add unnecessary factors to be considered, with the exception of reason 
(b) above which is a valid reason for the use of subroutines in most 
computer applications. 

By careful construction of the statement numbering scheme, available 
in languages such as FORTRAN, the programmer of the moderate to small size 
simulation can modularize his program without the additional consideration 
of designating common storage and the other difficulties experienced when 
programming subroutines. 



ADVANTAGES 

The above technique allows the designer to use the familiar computer 
languages, such as FORTRAN, in place of special simulation languages. 
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It also alleviates the requirement of either providing the reader with a 
description of a special language or causing him to go to another source 
to interpret the program. Using a "standard" language, such as FORTRAN, 
the designer can reasonably assume that the reader needs little or no 
explanation. 



FLEXIBILITY 

The use of modularization with statement numbers gives all the 
inherent flexibility observed in the special languages when the simulation 
is of such a magnitude as to allow compilation in one pass. When simu- 
lations are in the design and programming stages, the use of the above 
technique can allow as much flexibility as the simulation requires. Exit 
from the blocks when programming in FORTRAN can be accomplished at any 
logical point with nothing more than a simple GO TO or COMPUTED GO TO 
statement, and thus the logic flow is easily accomplished by this method. 
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4. SYSTEM VARIABLES, COORDINATES AND NOTATION 

Certain special simulation languages provide for dynamic allocation 
of storage for tables which facilitates the designation of data in a 
flexible and expandable form. This is not necessarily required by the 
moderate size simulation. The designer of these moderate size simulations 
may find the techniques as explained below more desirable, since his 
problem is not generally one of storage limitations. The criterion for 
the dynamic versus preset storage decision is felt to be that of program 
size. In programming the moderate size simulation the programmer may 
find that constructing a scheme for naming variables may create easier to 
work with variables than the complex tabular form of the special simu- 
lation language. This was found to be true in the programming of the 
example simulation, MULNUCl, of this thesis. 

The general approach to the naming of variables in the example 
program was that of using vectors to represent a given parameter, each 
component of the vector being representative of the value of that parameter 
with respect to the unit concerned. An example, that of the X-coordinate 
of the i t ^ destroyer, being DDX(I) . Once the scheme is understood the 
programming moves along without a great amount of thought required by the 
programmer as far as variable names are concerned. Certain variables 
must be set aside as dummy or temporary and these logically take on forms 
such as: ITEMP, TEMP3, DUMMY(I), and forms that immediately classify them 
in this category. 

In general, several different coordinate systems are required in 
simulations. In the war game an overall "area of play" must be established. 
This can be either rectangular or polar, two or three dimensional. The 
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axes of the rectangular coordinate system are not necessarily graduated 
in the same units. In the example program, for instance, it will be 
seen that the third axis (depth) is dimensioned in feet, while the two 
major axes are dimensioned in yards. Further, the third axis has the 
normally negative direction established as positive. These, perhaps 
unorthodox, measures are taken to facilitate programming; however, they 
must be spelled out in the documentation of the simulation to prevent 
possible misunderstanding. 

Many simulations require more than one coordinate system to be 
employed. The overall play is perhaps in a rectangular coordinate system, 
while range and bearing information may be required during the play of 
the game. This will necessitate the incorporation of an overlay of one, 
or perhaps several, polar coordinate systems upon the base system. If 
a "close-up" view is required during the play a translation and/or 
expansion to another rectangular system may be required. It can now be 
seen that, in general, several cooidinate systems will be used in a war 
game type of simulation. Therefore, a plan for the designation of these 
various coordinate systems and their respective transformations must be 
established early in the formulation of a simulation. 
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5. MULNUC1 - AN EXAMPLE 



MULNUCl is an on-line, real time simulation used as an example of 
the application of ideas presented in sections one through four. The 
simulation had its beginning at the Naval Radiological Defense Laboratory, 
Hunters Point, in the summer of 1965. During a six-week tour of the 
Naval Radiological Defense Laboratory, it was found that little had been 
done in the exploration of tactics and possible reactions of surface 
antisubmarine destroyers exposed to a self-inflicted multiple burst 
nuclear environment. At this time the simulation used as an example in 
this thesis had its beginning. In the example program, MULNUCl, all 
classified input parameters have been assigned fictitious values so that 
the computer program, as presented in this thesis, could remain unclassified. 

In 1965, Lieutenant J. E. Johnson programmed the on-line display, 

Display Data Corporation model DD 65, using a rather simple simulation 
situation [23. He did make a contribution in the form of an advancement 
in the techniques of on-line programming of simulations. Having observed 
a demonstration of Johnson’s program, it was felt that the technique of 
on-line display would be ideal for the envisioned program, MULNUCl. 

Several links were missing in the chain necessary to put the envisioned 
program on-line. The first link was a requirement for a routine to gen- 
erate circles of arbitrary size and location. This was accomplished with 
subroutine CIRCLE (see Appendix III). After completing subroutine 
CIRCLE, attention was turned to the necessity for a random number gen- 
erator capable of generating several types of random variables. The 
distributions of random variables required were uniform, normal, and 
circular normal. These generators were written in the form of subroutines 
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— apparent communications 

FIGURE 1 
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UNIFORM, NORMAL, and ERROR. The subroutine RANVAR is the basic random 
number generator called by these subroutines in the generation of their 
respective random variables. The above routines completed, only the 
communication routines required to link the Control Data Corporation 
1604 and 160 computers remained. This requirement was satisfied by sub- 
routines DCIRCLE, DTRACK, PARAMS and DSTATUS (see Appendix III and 
Acknowledgements) . 

At this point the preliminary work was complete and the formulation 
of the initializing, iterative loop, and critique portions of the simu- 
lation was begun. These three basic steps, as discussed in section two, 
were incorporated into an executive control block. 

EXECUTIVE CONTROL 

The Executive Control routine is flow charted in Appendix II and 
consists of three major parts. The first of these parts is the 
initializing portion. It is made up of four blocks: 

1. Inputs 

2. Set Constants 

3. Initialize 

4. Enter Input Changes 

In this subsection we shall consider the first three, leaving the 
latter for discussion in the subsection titled Man Machine Interface. 

The Inputs block is the one in which "standard" or nominal input 
parameters are set. These parameters are listed below. 

Number of Destroyers 
Time Factor 

Destroyer and Submarine Maximum Speed 
Submarine Hull Parameters 



18 



Initial Positions 



Initial Courses and Speeds 
Nominal Yield of Nuclear Weapons 
Depth of Thermocline 
Maximum Range of Weapons 
Detonation Parameters for Weapons 
Random Number Generator Initializer 

Any of these parameters may be changed in the enter input changes block 
(see Man Machine Interface subsection) . These parameters were chosen 
as the minimal requirements necessary to produce a simulation that has 
some realism and yet is not too complex. The structure of this program 
i6 such that any block can be expanded to include more parameters, 
thereby creating a more realistic simulation. 

Set Constants in a block used, as the name implies, to Initialize 
non-changeable inputs. In this block, all the logic indicators are set 
to orient the game. All damage and radiation levels are set at zero. 
The indices for tracking, firing, and sonar contact are set at zero. 
This is easily followed by cross referencing Appendices I and IV. 

The Initializing block begins by initializing the random number 
generator subroutine and then calculates the following parameters: 

Water Temperature Gradient 
Submarine Crush Depth 
Operational Depth of Submarine 
Wind Direction and Velocity 
Minimum Safe Range of Weapons [6] 

Effective Sonar Range 
Sea State 
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These calculations are straight forward and can easily be followed by 
cross referencing Appendices 1 and IV. The block also presents input 
parameters of interest to the destroyer team, submarine team (if selected) , 
and critique routine. 

Subsurface nuclear bursts are divided into four classifications: 

1. Very Shallow 

2. Shallow 

3. Deep 

4. Very Deep [4] 

The criteria for selection of classification are depth of burst and yield. 
The determination of classification of burst is made at this point in the 
program and an index IDEEP is set (see Appendix I) . The four matrices of 
output data are filled with negative zero, since negative zero is pro- 
grammed not to print on the display. Finally, if the role of the sub- 
marine is to be played by the computer, the basic strategy of the sub- 
marine is randomly determined (see Submarine Logic Model) . 

The next major part of the Executive Control Routine is the iterative 
loop. It is made up of eight blocks: 

1. Plot Positions 

2. Display Data 

3. Plot Generator 

4. Interactions 

5. Radiation Model 

6. Enter Changes 

7. Submarine Logic Model 

8. Time Loop 

The first three of these will be considered now, while the remaining blocks 
will be explained in later subsections. 
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The Plot Positions block is utilized to transmit output information 
to the left tube of the DD 65 (see Figure 2). This output information 
consists of the following data: 

Destroyer Tracks 
Sonar Contact Plots 

Destroyer Courses, Speeds and Coordinates 
Orientation and Size of Area Displayed 
Important Messages to the Player 

This information presentation is covered in more detail in the section 
on Man Machine Interface. 

The Display Data block performs the same function with respect to 
the right tube of the DD 65 (see Figure 3). Figure 3 lists the data 
displayed by this block and for this reason it will not be listed at 
this point. 

The final major part of the Executive Control Routine to be con- 
sidered at this time is the Critique block. This block critiques the 
simulation by performing several tasks. The entire simulation is recorded 
on magnetic tape and can be reviewed at a later time. If any significant 
action or interaction occurs game time and the nature of the action/ 
interaction are recorded by the Critique I block. A narrative print out 
is made at the conclusion of the simulation by the Critique II block. 

To correlate this information a graph plot of the tracks of destroyers, 
submarines, and all pools and clouds of radiation is made. Detail of 
these operations is documented in Appendix IV. 

PLOT GENERATOR 

The Plot Generator block advances all participants each time step 

(see Appendix II) . The block then determines if there are any clouds or 

pools of radiation. If there are clouds, they are advanced using wind 
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INPUT PAR AME I ERS 

NUMBER OF DESTROYERS ... 2 

SIZE OF ASROC WARHEAD ... 2.0 KT 

DEPTH OF THEKMOCLINE ...120.0 FT 

MAXIMUM RANGE OF ASROC . 9000.0 YARDS 

DEPTH OF BURST ' 700.0 FT 

SINK RATE OF WARHEAD ... 18.0 FT/SEC 

TIME FACTOR 5.0 
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FIGURE 3 



velocity and direction to determine motion, and the radius is computed 
[ l] • Pools are considered to be stationary since their drift is negligible* 

INTERACTIONS 

Interaction is the largest block and for clarity has been broken 
into sub-blocks. These sub-blocks are: 

1. Sonar Contact Model 

2. Contact Tracking Model 

3. Weapon Firing Model 

4. Evaluation Model 

The interactions block is constructed in the form of a loop that considers 
the interactions of each unit in succession. The organization is well 
documented in Appendices II and IV and will not be further explored at 
this time, however, the details of the four sub-blocks listed will be 
considered below. 



SONAR CONTACT MODEL 

The Sonar Contact Model uses a ray path theory detection scheme in 
a deterministic manner [3]. This deterministic detection range then has 
a variance superimposed upon it. The net result is a fairly realistic 
sonar detection model. The major limitation of this model is that it 
only handles the constant temperature gradient case. The inclusion of 
other gradients causes the sonar detection problem to assume a much more 
complex nature. Included in the model is a degradation of sonar range 
due to excessive destroyer speed. The range and bearing given as outputs 
from this model have range and bearing errors included. The data, with 
these errors, is then utilized by the tracking and firing models. Doppler 
is also calculated in the model and sent to the display as: 
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1. Up Doppler 

2. Down Doppler 

3. No Doppler 



CONTACT TRACKING MODEL 

The Contact Tracking Model uses as input data the output of the 
Sonar Contact Model. A simple criteria, requiring three consecutive 
marks from the sonar model, is used to distinguish non-contacts from 
contacts. Once three consecutive marks are received, the model determines 
the contact course and speed. This is done with a simple no parameter 
track model that considers at least three but no more than five marks. 

The course and speed of the contact are determined from the first and 
last of these marks. This track model is unsophisticated but the error 
induced in the output is fairly realistic. This model is of the first 
that should be improved upon if more work is to be done on this simulation. 
The output of this model is in the form of contact course and speed. This 
model provides dual routing, dependent upon the track. This will be 
discussed further in the subsection on Block Sequencing. Determination 
of whether or not the target is in firing range is made just before exiting 
the routine. 



WEAPON FIRING MODEL 

The Weapon Firing Model takes the last position of the contact from 
the Sonar Model, the contact course and speed from the Tracking Model, 
and then determines a firing solution. The range is considered, with 
time of flight and sink time, and the time of burst is determined. A 
dead reckoning position of the target is computed from the track data 
and this position becomes the aim point of the weapon. The model then 
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calls subroutine ERROR from which the true fall of shot is determined. 



The weapon is given a reliability check in the model and if this test is 
failed, the player will be notified that the weapon has misfired (see 
Man Machine Interface subsection) . 

EVALUATION MODEL 

The evaluation model initializes the pool and cloud of radiation 
created by the subsurface nuclear burst. To accomplish this the model 
takes data from the firing model for the location of ground zero and 
data from the inputs block for yield, depth of burst, and type of burst. 
The radius of the radioactive pool is then determined [4] (the cloud 
parameters are computed in the plot generator model) . 

A simple criteria for damage to the submarine is used. The lethal 
range is determined using submarine hull parameters, yield of warhead, 
and submarine depth as received from the submarine logic model. Slant 
range to the burst from the submarine is computed. If the submarine is 
within the lethal range, damage is 100%. If the submarine is outside a 
radius equal to twice the lethal range, damage is zero. Values of sub- 
marine damage between zero and 100% are computed by a linear relation- 
ship, then, if at any time the submarine's total damage reaches the 75% 
level, the game is terminated with the submarine considered as having 
been sunk. 



RADIATION MODEL 

The Radiation Model determines if any weapons have been detonated. 
If none have been detonated, the block is bypassed and the program con- 
tinues. If weapons have been detonated, the model computes the distance 
of each destroyer from all pools and clouds of radiation. The location 
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and size of all pools and clouds is received from the plot generator 
model along with the location of the destroyers. The model then determines 
if the destroyers are within the perimeter of any pool or cloud. If thi6 
condition exists, the radiation level in the cloud or pool is calculated 
[1> 4]. The total radiation being received by each destroyer is then 
calculated. The radiation rate and total radiation dose for each 
destroyer are computed and sent to the display (see Man Machine Interface 
subsection) . 



SUBMARINE LOGIC MODEL 

The submarine can be controlled in two ways: 

1. By a submarine team. 

2. By the computer. 

This decision is made in the enter input changes routine. The program 
is such that the computer will play the role of the submarine unless the 
variable ISUB is set equal to one by the enter input changes routine. 

If ISUB is set equal to one, control of submarine depth, course, speed, 
and weapon firing is turned over to the submarine team. This team will 
receive passive sonar bearings and screw beat information from the 
console typewriter of the CDC 1604. They will be able to control the 
movements of the submarine and its weapon firing by means of the console 
typewriter and selective jump keys. The program will query selective 
jump key number two, once each cycle, to determine if orders to the 
submarine are to be received. If selective jump key two is set the 
computer will request orders (see Man Machine Interface subsection) . 

The weapons (torpedoes) can be fired by setting selective jump key three 
on the CDC 1604 console (see Man Machine Interface subsection) . Under 
this condition the simulation becomes a conflict between two teams: 
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1. A submarine team using the CDC 1604 console. 

2. A destroyer team using the DD 65 display console. 

If ISUB is unchanged by enter input changes (inputs block sets 
this variable equal to zero) , the moves of the submarine are controlled 
by the computer. This being the case, two basic initial tactical sit- 
uations are available to the player. The first of these places the 
submarine on the surface, at the origin of the playing area. The sub- 
marine knows that he has been sighted and the game proceeds. The second 
situation has the submarine randomly located in the upper half of the 
playing area. In this case the submarine's position is not known by the 
player. It should be noted that this is the initial situation if the 
submarine is to be controlled by a submarine team. 

Initial situation one is selected by setting INITIAL equal to zero, 
while situation two is selected by setting this variable equal to one. 
Having selected the initial situation the computer then selects one 
of three basic strategies: 

1. The submarine runs for it. 

2. The submarine tries to transit between the two destroyers. 

3. The submarine tries an end run, flanking the two destroyers. 
These strategies are easily followed in Appendix II and will not be 
explored further at this point. 

SEQUENCING OF PROGRAM BLOCKS 

The blocks in this simulation are of two basic types: 

1. A single point of exit . 

2. Multiple points of exit. 

The block that has a single exit point might be called a standard block, 
in that, the block is called upon to perform a computation but no 
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branching of logic is done within the block. The multiple exit block 
is one in which the program is routed differently depending upon logic 
decisions made within the block. The sonar contact block is an example 
of this type. In this block the routing depends upon the results of a 
sonar search. If contact is made, the block exists to the contact 
tracking block. If no contact is made, the block exits to consider the 
next destroyer. 

The executive control routine flow chart in Appendix II illustrates 
the time sequencing of the major program blocks. The simulation is 
delayed at four points in the program which are: 

1. Enter input changes routine. 

2. Enter changes. 

3. Submarine logic block. 

4. The time loop. 

The first of these interruptions takes place only during the initializing 
portion. At this point any change to the input parameters is made. The 
second of these interruptions is made once during each time step. This 
is the point at which destroyer team changes are sent to the CDC 1604. 

The third of these interruptions takes place if the situation using the 
submarine team has been selected. In this case the simulation may be 
interrupted every time step to allow the submarine team to enter changes. 
The fourth of these interruptions occurs each cycle and maintains the 
time stepping interval. 

It will be noted that in the executive control routine, each block 
is considered in turn, no block is bypassed. In the interactions block, 
it will be noted that, sub-blocks are not always considered. No sonar 
contact by the sonar contact block causes the contact tracking block to 
be bypassed. The same is true of the weapon firing block if the tracking 
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block does not generate a satisfactory track. 

MAN MACHINE INTERFACE 

This subsection is concerned with communications, both into and 
out of the computer. These man-machine communications fall into four 
types: 

1. Data to and from the destroyer team. 

2. Data to and from the submarine team. 

3. Data to the various modes of the critique routine. 

4. Input changes. 

The first type breaks into three pares: 

1. Right tube information. 

2. Left tube information. 

3. Changes sent to the CDC 1604. 

The right tube gives the destroyer team data in tabular form as illustrated 
in Figure 3. The left tube will display the tracks of the destroyers, 
any sonar contacts, and all pools and clouds of radiation. Also included 
in the display on the left tube is a series of windows in which data can 
be displayed. The windows are numbered as shown in Figure 2. Window 
data assignments are as listed below. 



Window 


Data Assignment 


Units 


1 


X-coordinate of destroyer 1 


100 yards 


2 


X-coordinate of destroyer 2 


100 yards 


3 


Y-coordinate of destroyer 1 


100 yards 


4 


Y-coordinate of destroyer 2 


100 yards 


5 


Course of destroyer 1 


Degrees true 


6 


Speed of destroyer 1 


Knots 
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Window 


Date Assignment 


Units 


7 


Course of destroyer 2 


Degrees true 


8 


Speed of destroyer 2 


Knots 


9 


X-coordinate of left tube center 


100 yards 


10 


Y-coordinate of left tube center 


100 yards 


11-14 


Available for flash messages 


Alfa-numeric 


15 


Not used 




16 


Radius of display on left tube 


1000 yards 



Windows 5-10 and 16 are controllable from the DD 65 console by means of 
a discrete digital type control system. In this manner the player is able 
to change the destroyer course and speed, or "zoom" the display in on any 
point in the area of play. Let us first consider windows five through 

eight. 

On the DD 65 console (see Figure 4) there are buttons labelled 
CONN DD 1 and CONN DD 2. By depressing one or both of these buttons the 
player is given control of the course and speed of the destroyer or 
destroyers selected. Near the button just selected is a group of four 
buttons (see Figure 4) labelled RIGHT, LEFT, FAST/UP, and SLOW/DOWN. 
Depressing one of these buttons will cause the appropriate variable in 
windows five through eight to change in the desired direction. As an 
example, if the player depresses both CONN buttons and then the button 
labelled RIGHT - both destroyers will commence a turn to the right. 

They will continue this turn until the button labelled RIGHT is released. 
The player can "come to a course" by depressing the correct button until 
the desired course is displayed in the respective window. The same 
procedure is used for changing speed. It should be noted that the 
windows can be changed to values that are unacceptable, such as -5 knots, 
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in this case the program will cause this value to be changed back into 
the acceptable range during the following cycle. The value of course 
can be increased to values greater than 360 degrees, in which case the 
program will correct to the acceptable value. For example, if the player 
increases the course to 390 degrees, the program will convert this to 
030 degrees during the next cycle. 

The same basic procedure is utilized to "zoom" the display to any 
location desired. First the player selects the button labelled SHIFT 
(see Figure 4). He now has control of windows nine and ten. Depressing 
the button labelled RIGHT will cause the display to shift to the right, 
LEFT accomplishes the same action but to the left, and similarly with UP 
and DOWN. To zoom in, the rotary switch in the upper right hand corner 
of the console is used (see Figure 4). This switch has six positions 
labelled 4, 8, 16, 32, 64, and 128. Selection of one of the six positions 
will cause the radius of the displayed area to be that of the value 
selected, in thousands of yards. As an example, selecting SHIFT and 
changing windows nine and ten to 240 and 150 respectively, then setting 
the rotary switch to eight, will cause the area centered at (24000,15000) 
with a radius of 8000 yards to be displayed on the left tube. Note that 
window 16 will show the value eight, while windows nine and ten will show 
240 and 150 respectively. 

Windows 11 through 14 are used to send the player the following 
flash messages. 

SONAR CONTACT 
ASROC FIRED 
ASROC MISFIRE 
SUB SUNK 
DD SUNK 
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I- 








a 

fa 
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TARGET IN RANGE 



TARGET TOO CLOSE TO SHOOT 
TARGET OUT OF RANGE 

Another method of communications is available to the destroyer team. 

By depressing the DD 1 FIRE button, destroyer number one fires an ASROC. 
The same procedure is used with the DD 2 FIRE button. 

The next type of communications available is that of the submarine 
team. This is, of course, non-existent if the option is chosen in which 
the submarine is played by the computer. The submarine team will receive 
messages each cycle containing bearing and screw beat information. The 
submarine team may then choose to maneuver the submarine by setting 
selective jump key number two on the console of CDC 1604. The computer 
will then type COURSE ORDERS on the console typewriter. This is the 
indication that the computer is ready to receive course changes. If no 
change is desired, the old course is typed in. If a change is requested, 
the new course is typed in. The course typed in should be of the form 
090. followed by a carriage return. This will change the course and the 
computer will return NEW SUB COURSE 090. This completes the course change 
cycle and the computer will then type SPEED ORDERS, the same procedure 
is used to enter speed changes. Upon completion of the speed entry the 
computer will return NEW SUB SPEED 15, followed by DEPTH ORDERS. The new 
depth is now entered, and the computer will return NEW SUB DEPTH 1050. 

The routine is now finished and the program continues. 

The only other action available to the submarine team is that of 
firing torpedoes. This is accomplished by depressing selective jump key 
number three until the typewriter returns TORPEDO FIRED. At this time 
the selective jump key should be returned to the normal position, unless 
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another torpedo is desired. If a hit is scored, the typewriter will 
return DESTROYER SUNK. When both destroyers are sunk the typewriter 
will return GAME OVER. At the beginning of the game, the typewriter 
will give the submarine team the following information: 

1. Maximum speed available to submarine. 

2. Maximum depth allowable. 

The third type of communication is with the various critique routines. 
Critique is accomplished in three ways: 

1. A recording of all information on the DD 65 display is recorded 
by the tape unit near the DD 65. 

2. A graph of the tracks of the DD's, submarine, and all radiation 
is made on tape unit eight of the CDC 1604 (the graph has game time 
recorded by each mark to aid in correlating with the various other parts 
of the critique routine). 

3. A critique of all important events and their time is recorded 

on tape unit five of the CDC 1604 for print out at the conclusions of the 
game. 

These three methods of critique, if correlated, will give an excellent 
"replay" of the simulation. Any communication with the program other than 
listed above will be accomplished as described by Leach and Perrella [5]. 

TIMING 

The timing of this simulation is done by means of a time loop block. 

In this block, the contents of cell 5006B in the CDC 1604- is tested and 
stored as ICLOCK. A variable, NEXT, is generated as the sum of ICLOCK 
and ISTEP, the stepping interval. ISTEP is determined from another variable 
TFACTOR that is equal to one for real time. TFACTOR equal to three \;ould, 
for example, cause the game to run at three times real time. When ICLOCK 
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becomes greater than or equal to NEXT the loop is exited and the simu- 
lation continues. This procedure is easily followed in Appendices 
II and IV. 
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6. CONCLUSIONS AND ACKNOWLEDGEMENTS 



The considerations made in the formulation and construction of 
this simulation have made some observations possible. Only the basic 
modular structure of the simulation has been considered in detail. Each 
individual modular block has been designed in as simple a form as 
possible while maintaining some degree of realism. The simulation is 
sound in its general organization. New program blocks may be substituted 
making the simulation as realistic as desired. It is hoped that this 
simulation will be played with more sophisticated models and on-line 
equipment of greater capacity so that doctrine and tactics in the area of 
self-inflicted nuclear environment may be explored. 

It was found that the general purpose computer language was com- 
pletely satisfactory for the construction of this simulation. The mod- 
ularization technique made the logical organization of the simulation 
straight forward and is recommended for use in future simulations. 

The author wishes to express his appreciation to Professor 
Mitchell L. Cotton and Professor Alvin F. Andrus for their aid and 
encouragement in the preparation of this thesis. In addition, the author 
would like to thank Miss Patricia Hoang for her assistance in programming 
the linking subroutines and LCDR Richard E. DeWinter for his comments 
in proofreading the manuscript. 
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APPENDIX I 



LIST OF VARIABLES 

This Appendix contains a listing of the variables used in the 
simulation MULNUCl, arranged in alphabetical order. 



Variable 


Definition 


AROCMAX 


Maximum range of the ASROC, in yards. This is an 
input parameter, set equal to 9,000 yards by the inputs 
block and can be changed in the change inputs block. 


AROCMIN 


Minimum safe range for the ASROC, in yards. This value 
is computed in the initializing block, and is a function 
of warhead size. 


B 


Last bearing the submarine held of the nearest destroyer, 
in degrees true. 


BNPTS 


A dummy variable used in the contact tracking block 
for the determination of contact speed. 


CENTERB 


The course such that the submarine will split the 
channel between the destroyers, in degrees true# 


CLOUDR(I) 


t'Vi 

Radius of the i n cloud of radiation, in yards. 


CLOUDX(I) 


X-coordinate of the i fc h cloud of radiation with respect 
to the main coordinate system, in yards. 


CLOUDY (I) 


Y-coordinate of the i fc ^ cloud of radiation with respect 
to the main coordinate system, in yards. 


CONTB(I) 


Bearing of sonar contact of i*"* 1 destroyer, in degrees 
true. 


CONTC(I) 


Contact course, as computed by the contact tracking 
block, in degrees true. 


CONTR(I) 


Sonar range to contact, as measured by the i* 1 ' 1 destroyer. 


CONTS(I) 


Contact speed as computed by the contact tracking model, 
in knots. 


CONTX(I) 


X-coordinate of the i^ 1 destroyer's contact, in yards, 
with respect to the main coordinate system. 


CONTY(I) 


Y-coordinate of the i c ^ destroyer's contact, in yards, 
with respect to the main coordinate system. 
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Variable 


Definition 


D 


Difference between the true bearing to the submarine 
from the destroyer and the submarine's true course, 
in degrees. This is used in the determination of 
doppler. 


DA 


Absolute value of D. 


DAMAGE 


Percent damage to the submarine from the current 
detonation. 


DAMAGET 


Cumulative damage to the submarine, in percent. 


DDC(I) 


Course of the i^ destroyer, in degrees true. 


DDS(l) 


Speed of the i tn destroyer, in knots. 


DDSMAX(I) 


Maximum speed available to the i^ 1 destroyer, in knots. 


DDX(l) 


X-coordinate of the i^ destroyer, in yards, with 
respect to the main coordinate system. 


DDY(I) 


Y-coordinate of the i^ 1 destroyer, in yards, with 
respect to the main coordinate system. 


DETR 


Detection range, in yards. This is a random variable 
with mean DETRM and normally distributed with sigma 
of .3 times DETRM. 


DETRM 


Mean detection range, in yards. A function of GRAD 
SUBD. 


DISTC(I,J) 


Distance of the i t destroyer from the center of the 
jth cloud of radiation. 


DISTP(I,J) 


Distance of the x destroyer from the center of the 

jth pool of radiation. 


DOB 


Depth of burst, in feet, of ASROC warhead. 


DR 


Advance of the destroyer considered, in yards, between 
marks as considered by the tracking model. 


DTIME 


Total time delay from the firing of an ASROC and the 
detonation, in seconds. This includes time of flight 
and time of sinking. 


DUMMY 


A dummy variable used throughout the program. 


DX 


East-west advance of the destroyer considered, in yards, 
between marks as considered by the tracking model. 
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Variable 


Definition 


DY 


North-south advance of the destroyer considered, in 
yards, between marks as considered by the tracking 
model. 


GSR 


Effective sonar range, in yards. Computed by the 
initializing block. 


GRAD 


Water temperature gradient, in degrees per hundred 
feet of depth. 


GTIME 


Game time, in minutes. At the start of the game 
GTIME is zero. 


GZX 


X-coordinate of ground zero for the detonation 
considered, in yards, with respect to the main 
coordinate system. 


GZY 


Y-coordinate of ground zero for the detonation 
considered, in yards, with respect to the main 
coordinate system. 


HULL 


Submarine hull thickness, in inches of steel. 


IAS ROC 


Hull number of the destroyer that fired the last 
ASROC. If no ASROC has been fired, IASROC is zero. 


ICLASS(I) 


Sonar classification of the i^ destroyer's contact. 
0 - no contact, 1 - possible submarine, 2 - probable 
submarine. 


ICLOCK 


Contents of location 5006B in the CDC 1604. The CDC 
1604 steps this cell once every second. 


ICIRCX(I) 


A dummy vector of X-coordinates of points to transmit 
circles to the display. 


ICIRCY(I) 


A dummy vector of Y-coordinates of points to transmit 
circles to the display. 


ICONTB(I) 


Fixed point version of CONTB(I) , used to transmit to 
the display. 


ICONTC(I) 


Fixed point version of CONTC(I) , used to transmit to 
the display. 


ICRIT# 


A series of critique indicators. 0 - pass. 


1 


1 - sonar contact 


2 


1 - ASROC fired 
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Variable 


Definition 



ICRIT# 



3 


1 - torpedo fired 


4 


1 - ASROC misfire 


5 


1 - target in range 


6 


1 - sub sunk 


7 


1 - destroyer sunk 


8 


1 - target too close to shoot at 


9 


1 - target out of ASROC range 


11 


1 - game is a draw, submarine escaped 


12 


1 - submarine wins, destroyer 1 sunk with submarine 
escaping 


13 


1 - submarine wins, destroyer 2 sunk with submarine 
escaping 


14 


1 - submarine wins, both destroyers sunk 


15 


1 - destroyers win, submarine sunk by destroyer 1 


16 


1 - destroyers win, submarine sunk by destroyer 2 


17 


1 - submarine wins by transiting between the destroyers 


IDDC(I) 


Fixed point version of DDC(I), used to transmit to 
the display* 


IDDS(I) 


Fixed point version of DDS(I) , used to transmit to the 
display. 


IDDX(I) 


Fixed point version of DDX(I) , used to transmit to the 
display. 


IDDY(I) 


Fixed point version of DDY(I), used to transmit to the 
display. 


IDEEP 


An index used to indicate the classification of nuclear 
burst. 1 - very shallow, 2 - shallow, 3 - deep, 4 - 
very deep. 


IDOPLER(I) 


An index used to indicate the doppler of the i*"* 1 
destroyer's contact. 0 - no doppler, 1 - up doppler, 
2 - down doppler. 
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Variable 



Definition 



I END 
IFIRE 

ILOGIC 

IONE 

INITIAL 

IR 

IRANDOM 

ISHOOTl 

ISHOOT2 

ISOL(I) 

ISS 

ISTEP 

ISTRAT 

ISUBC 

ITURN 

1X0 

IXDD(I,J) 



An index used to indicate game over. 0 - game not over, 

1 - game over. 

An index used to indicate if the submarine has fired a 
torpedo. 0 - torpedo not active, 1 - torpedo still 
active. 

A logic index used in the submarine model. 

The lowest destroyer hull number. Equal to 1 at the 
beginning of the game. If destroyer number one is sunk 
the IONE is equal to two. 

An index of the initial situation. 0 - submarine at 
the origin, 1 - submarine randomly distributed in 
the upper half of the playing area. 

Radius of display area, in thousands of yards. 

A random number selected to initialize the random 
number generator. This number must be an odd integer 
in the interval 1 to 67108863. 

A flag used to activate the firing sequence for 
destroyer 1. 0 - do not shoot, 1 - shoot. 

A flag used to activate the firing sequence for 
destroyer 2. 0 - do not shoot, 2 - shoot. 

A progressive index of the quality of the firing 
solution the i** 1 destroyer has on its target. 0 - no 
solution, thru 5 - best solution. 

Sea state. 

Fixed point version of TSTEP. 

Submarine basic strategy when the submarine is controlled 
by the computer. 0 - run for it, 1 - go up the middle, 

2 - end run. 

Fixed point version of SUBC. 

Logical index that records the submarine's initial 
turn. 0 - left turn, 1 - right turn. 

X-coordinate of center of displayed area, in hundreds 
of yards with respect to the main coordinate system. 

A dynamic table of track data, recording the X-coordinate 
of the j*-* 1 mark of the ith destroyer. The maximum value 
of j is eight. Entries are in yards, with respect to the 
main coordinate system. 
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Variable 


Definition 


IXDDl(l) 


This variable is equal to IXDD(1, I) . 


IXDD2(I) 


This variable is equal to IXDD(2, I). 


IXSUB(I.J) 


A dynamic table of track data, recording the X- 
coordinate of the jth mark of the i c h destroyer's 
sonar contact. The maximum value of j is eight. 

Entries are in yards, with respect to the main coordinate 
system. 


IXSUBl(I) 


This variable is equal to IXSUB(1,I). 


IXSUB2(I) 


This variable is equal to IXSUB(2,I). 


IYO 


Y-coordinate of center of displayed area, in hundreds 
of yards with respect to the main coordinate system. 


IYDD ( I , J) 


A dynamic table of track data, recording the Y-coordinate 
of the jth mark of the ith destroyer. The maximum 
value of j is eight. Entries are in yards, with respect 
to the main coordinate system. 


IYDDl(I) 


This variable is equal to IYDD(1»I) . 


IYDD 2 (I) 


This variable is equal to IYDD(2,I). 


IYSUB(I.J) 


A dynamic table of track data, recording the Y-coordinate 
of the mark of the i fc ^ destroyer's sonar contact. 

The maximum value of j is eight. Entries are in yards, 
with respect to the main coordinate system. 


IYSUBl(I) 


This variable is equal to IYSUB(1, I). 


IYSUB2(I) 


This variable is equal to IYSUB(2, I). 


IWINDD 


Fixed point version of WINDD, used to transmit data to 
the display. 


MARKS (I) 


f*h 

Number of continuous marks, up to five, the i 
destroyer has on his sonar contact. 


N 


A dummy variable. 


NCONBER 


Number of constant bearing the submarine has on the 
nearest destroyer. 


NEXT 


Time of the next cycle, in seconds of computer time. 


NMARKS(I) 


Same as MARKS(I) except NMARKS(I) does not stop at five. 
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Variable 


Definition 


NOSHOOT 


Index to limit destroyers to one active weapon at a 
time. 0 - there are no active weapons, alright to 
shoot, 1 - there is an active weapon, cannot shoot. 


NPTS 


Index used to control display data until there are 
eight point available on destroyer tracks. 


NPTSl 


A dummy variable. 


NPTS 2 


A dummy variable. 


NRDD 


The highest destroyer hull number. Equal to 2 at the 
beginning of the game. If destroyer number two is 
sunk then NRDD is equal to one. 


NSHOTS 


Total number of ASROC' s fired during the game. 


PHIT 


A random variable used to determine if a torpedo, fired 
by the submarine, hit the destroyer. 


POOLR(I) 


Radius of the i^ 1 pool of radiation, in yards. 


POOLX(I) 


X-coordinate of the i^ pool of radiation, in yards, 
with respect to the main coordinate system. 


POOLY(I) 


Y-coordinate of the i*^ 1 pool of radiation, in yards, 
with respect to the main coordinate system. 


R 


Radius of the displayed area, in yards. 


REL 


A random variable used to determine ASROC reliability. 


RADDOSE(I) 


Total radiation dose the i 1 " destroyer has been exposed 
to during the game, in roentgens. 


RADRATE(I) 


The rate at which the i destroyer is receiving 

radiation, in roentgens/hr. 


RANDOM 


The random variable used to link the various random 
generators. 


RANGE 


A dummy variable used as a temporary storage in range 
calculations. 


RLETHAL 


The lethal range of the ASROC warhead, in yards. 


SAFETY 


Safety factor used to compute submarine maximum 
operating depth as a function of the crush depth of 
the hull. 


SB (I) 


Screw beat count of the i^ destroyer as measured by the 
submarine. 
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Variable 


Definition 


SIGMA 


A dummy variable used in the calling of normally 
distributed random variables. 


SR 


Sink rate of the ASROC warhead, in feet per second. 


SSB(I) 


t*h 

The actual beating of the submarine from the i 
destroyer, in degrees true* 


STRESS 


Yield strength of the steel used in the submarine 
hull, measured in thousands of pounds per square inch. 


SUBC 


The true submarine course, in degrees true* 


SUBD 


The actual submarine depth in feet* 


SUBDMAX 


The maximum allowable operational depth of the submarine, 
in feet. 


SUBS 


The actual speed of the submarine, in knots. 


SUBSMAX 


The maximum speed available to the submarine, in knots. 


SUBX 


The actual submarine X-coordinate with respect to the 
main coordinate system, in yards. 


SUBY 


The actual submarine Y-coordinate with respect to the 
main coordinate system, in yards. 


SUM 


A dummy variable used in the computation of radiation 
dose. 


TBURST (I) 


Time of detonation of the i^ 1 ASROC warhead, measured 
in game time. 


TEMP 


A dummy variable. 


TEMPI 


A dummy variable. 


TEMP2 


A dummy variable. 


TEMP3 


A dummy variable. 


TEMP4 


A dummy variable. 


TEMP5 


A dummy variable. 


TEMPCR(I,J) 


The radiation rate received by the i fc ^ destroyer from 
the j* 1 ' 1 cloud of radiation, in roentgens per hour. 


TEMPPR(I,J) 


The radiation rate received by the i c ^ destroyer from 
the cloud of radiation, in roentgens per hour. 
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Variable 


Definition 


TF ACTOR 


A time factor, 1 - real time, 2 - double time, .5 - 
half time. 


THERMO 


Depth of the thermocline, in feet* 


THETA 


A dummy variable used in the computation of courses and 
bearings. 


T INTER 


Time to intercept of a torpedo fired by the submarine, 
in minutes. 


TLOGIC 


A time storage point used to control timed logic. 


TOB 


Time of burst, this includes time of flight and sink 
time of the ASROC warhead, measured in minutes of game 
time. 


TOF 


Time of flight of the ASROC warhead, in seconds. 


TOS 


Sink time of the ASROC warhead, in seconds. 


TSTEP 


The time step for each cycle of the game, in seconds. 


TXDD(I,J) 


i.L. 

A temporary storage of the j past position X-coordinate 
of the i^ destroyer used in preparing track data for the 
display, measured in yards. 


TYDD(I, J) 


A temporary storage of the past position Y-coordinate 
of the it* 1 destroyer used in preparing track data for the 
display, measured in yards. 


TXSUB(I,J) 


A temporary storage of the j*** 1 past position X-coordinate 
of the ith destroyer's sonar contact used in preparing 
track data for the display, measured in yards. 


TYSUB(I, J) 


A temporary storage of the j t ^ 1 past position Y-coordinate 
of the i fc k destroyer's sonar contact used in preparing 
track data for the display, measured in yards. 


VEL 


The relative velocity of the torpedo fired by the sub- 
marine, with respect to the destroyer target, in knots. 


Wl 


A dunmy variable used to determine the destroyer that 
is closest to the submarine. 


W2 


A dummy variable used to determine the destroyer that 
is closest to the submarine. 


WINDD 


Wind direction in degrees true. 


WINDV 


Wind velocity, in knots. 
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Variable 


Definition 


X(I,J) 


The X-coordinate of the i tl1 destroyer's j*"* 1 continuous 
mark on its sonar contact, in yards, with respect to 
the main coordinate system. 


XDD(I, J) 


The X-coordinate of the i fc ^ destroyer j-1 steps ago, 
in yards, with respect to the main coordinate system. 


xo 


The X-coordinate of the center of the displayed area, 
in yards, with respect to the main coordinate system. 


XTEMP 


A dummy variable used for temporary storage of X-coordinates 


Y(I,J) 


The Y-coordinate of the i* 1 * 1 destroyer's continuous 

mark on its sonar contact, in yards, with respect to the 
main coordinate system. 


YDD(I, J) 


The Y-coordinate of the i tn destroyer j-1 steps ago, 
in yards, with respect to the main coordinate system. 


YIELD 


The yield of the ASROC warhead in kilotons. 


YO 


The Y-coordinate of the center of the displayed area, 
in yards, with respect to the main coordinate system. 


ZZX 


A vector of destroyer and submarine X-coordinates used 
by the graph plot routine in Critique I. 


ZZY 


A vector of destroyer and submarine Y-coordinates used 
by the graph plot routine in Critique I. 
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APPENDIX II 



LOGIC FLOW DIAGRAMS 

This Appendix contains a series of logic flow diagrams as listed 
below. 

Executive Control 
Plot Generator 
Interactions 
Radiation Model 
Submarine Logic Model 
Run for It 
Up the Middle 
End Run 

Submarine Team Control 
Sonar Contact Model 
Contact Tracking Model 
Weapon Firing Model 
Evaluation Model 
Critique I 
Critique II 

It will be noted that the statement number from the program listing is 
shown in the upper right hand corner of each symbol in the block diagram. 
This should aid in correlating this Appendix, with Appendix IV. 
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FLOW CHART SYMBOLS 



C 





A connector or terminal. 



An offpage connector* 




A predefined process or module/subroutine. 

A more detailed flow chart of this subroutine is al 
included. 





Input/output other than display. 



Decision. 



Processing, annotation. 



Display 
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APPENDIX III 



SUBROUTINES AND CDC 160 EXECUTIVE ROUTINE 

This Appendix contains as explanation of the subroutines used in 
the main program and the CDC 16G executive routine used to connect the 
CDC 1604 with the DD 65 display. A listing of these subroutines can be 
found in Appendix IV. 



SUBROUTINE RANVAR 

Subroutine RANVAR is used as the basic random number generator. It 
generates random floating point numbers in the interval zero to one that 
are distributed uniformly in that interval. The random variables used 
in the main program are called from either UNIFORM, NORMAL or ERROR 
subroutines which in turn call RANVAR for input. The generator is a 
simple fixed point division utilizing the remainder from the Q register 
as the random number. This number is then mapped into the zero to one 
interval. Only one input is required to initialize this generator, 
namely IRANDOM. 



SUBROUTINE UNIFORM 

Subroutine UNIFORM is a three argument subroutine used to generate 
uniformly distributed random numbers in any interval. The arguments are: 

1. The center of the interval. 

2. The half width of the interval. 

3. The output random number. 

An example would be UNIFORM (5.0, 2.0, SUBS) • This call would yield the 
variable SUBS, submarine speed, uniformly distributed in the interval 
with center at five knots and plus or minus two knots. 
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SUBROUTINE NORMAL 



Subroutine NORMAL is a three argument subroutine used to generate 
normally distributed random numbers with any mean and standard deviation. 
This normal distribution is generated by means of the sum of identically 
distributed (uniform) random variables. Twelve uniform numbers are used 
because: 

1. The truncation is not significant. 

2. Twelve reduces the formula to a summation and no division is 
required. 

The arguments to this subroutine are: 

1. The mean of the distribution. 

2. The sigma of the distribution. 

3. The output random number. 

SUBROUTINE ERROR 

Subroutine ERROR is a five argument subroutine used to generate 
circular normal distributed random variables. This subroutine is used 
to determine the true fall of shot given the aiming point and CEP or 
sigma of the distribution of fall or shot. The arguments are: 

1. The X-coordinate of the aim point. 

2. The Y-coordinate of the aim point. 

3. The sigoa of the circular normal distribution. 

4. The X-coordinate of the fall of shot. 

5. The Y-coordinate of the fall of shot. 

SUBROUTINE CIRCLE 

Subroutine CIRCLE is used to generate a series of points, 24 in 
number, every 15 degrees around the perimeter of a circle of predetermined 
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center and radius. The circle is then used in the formation of the 
circular clouds and pools of radiation displayed to the destroyer team. 
The arguments are: 

1. The X-coordinate of the center of the circle with respect to 
the main coordinate syste. 

2. The Y-coordinate of the center of the circle with respect to 
the main coordinate system. 

3. The radius of the circle with respect to the main coordinate 
system. 

4. The X-coordinate of the center of the display with respect to 
the main coordinate system. 

5. The Y-coordinate of the center of the display with respect to 
the main coordinate system. 

6. The radius of the displayed area. 

7. A vector of X-coordinates of the 24 points in the circle with 
respect to the display coordinate system. 

8. A vector of Y-coordinates of the 24 points in the circle with 
respect to the display coordinate system. 

SUBROUTINE DCIRCLE 

DCIRCLE is a subroutine used to transmit the circle coordinates, 
generated in subroutine CIRCLE, to the CDC 160 from the CDC 1604. The 
arguments to this subroutine are: 

1. ITRKNO - the track of circle number by which this particular 
circle can be designated. 

2. CHAR - a single letter or number in hollerith form that is to 
be displayed as one of the points in the circle 

3. In the program MULNUCl the letter c and p are used to dis- 
tinguish between clouds and pools of radiation. 
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4. NUMPTS - the number of points in the circle. 

5. IX - a vector of X-coordinates for the circle. 

6. IY - a vector of Y-coordinates for the circle. 

SUBROUTINE DTRACK 

Subroutine DTRACK is used to transmit track data to the CDC 160 
from the CDC 1604. This subroutine will send tracks of up to eight points 
to the CDC 160. The arguments to this subroutine are the same as those 
in DCIRCLE, in fact the two routines are identical with the exception of 
the allowable number of points. 

SUBROUTINE DSTATUS 

DSTATUS is a subroutine used to transmit information to the windows 
described in the section on Man Machine Interface. This information can 
be in two forms: 

1. Program variables in fixed point form. 

2. Messages in hollerith form. 

The arguments are: 

1. ITYPE - zero represents a numerical program variable is to be 
sent, while one indicates that an eight hollerith character word is to be 
sent. 

2. NWIND - the number of the window to which the variable or word 
is to be sent. 

3. IW - the field width of the variable. 

4. INAME - the name of the variable to be transmitted. 

5. IX - the X-coordinate of the lower left hand corner of the 
window. 

6. IY - the Y-coordinate of the lower left hand corner of the 
window. 
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It should be noted that windows are 128 display coordinates units long* 

If a message in the form of words is to be sent to the display and it is 
longer than eight letters (the length of one window) it can be sent by 
means of more than one window. These windows should be displaced by 128 
units in the X direction, thus the windows may form a continuous word of 
more than eight characters. 

SUBROUTINE PARAMS 

PARAMS is an eight argument subroutine used to query the CDC 160 
as to the contents of eight selected windows# This subroutine allows 
the main program, in the CDC 1604, to enter changes that have been made 
to the windows of the display by the player. This is the only method the 
player has of communicating with the program without interrupting the 
play. The eight arguments to this subroutine are: 
lA - the contents of window 5 
2A - the contents of window 6 
3A - the contents of window 7 
4A - the contents of window 8 
5A - the contents of window 9 
6A - the contents of window 9 
7A - the contents of window 16 

8A - the numeric value of a location in the CDC 160 that is 
controlled by the DDl FIRE and DD2 FIRE buttons. 



73 



APPENDIX IV 



PROGRAM LISTING 

This Appendix contains the computer program listing of the simu- 
lation MULNUCl. The program blocks are in numerical order and the logic 
can be followed by cross referencing this Appendix with Appendix I. In 
this program all classified input parameters have been assigned fictitious 
values so that the program, as presented in this thesis, could remain 
unclassified. 
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PROGRAM MULNUC1 MUL00020 
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